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Abstract 

We study large-scale structure formation in the presence of a quintessence component with zero speed 
of sound in the framework of Eulerian Perturbation Theory. Due to the absence of pressure gradients, 
quintessence and dark matter are comoving and can be studied as a unique fluid in terms of the total 
energy density contrast and the common velocity. In this description the clustering of quintessence 
enhances the linear term proportional to the velocity divergence in the continuity equation by a factor 
(1 + w)Q.Q/£l m . This is responsible for a rapid evolution of the growth rate at low redshifts, and 
modifies the standard relation between the velocity divergence and the growth factor. For the total 
fluid, the solutions for the linear growth function and growth rate can be written in integral forms 
and admit simple fitting formulae, as in the ACDM case. At second order in perturbation theory, we 
derive an explicit expression for the kernels F2 and G2. They receive modifications of the order of 
the ratio between quintessence and total energy density perturbations, which affect the correspond- 
ing tree-level bispectra. We finally compute the cumulative signal-to-noise in the power spectrum, 
bispectrum and reduced bispectrum, expected for departures from a ACDM cosmology both in the 
clustering and smooth quintessence scenarios. The reduced bispectrum, in particular, receives sensi- 
ble modifications only in the clustering case and can potentially be used to detect or rule out the model. 
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1 Introduction 



The simplest and most popular model of dynamical dark energy is quintessence, a single scalar field 
whose vacuum energy dominates the Universe driving its acceleration. Quintessence energy density 
varies with time and a way to distinguish it against a cosmological constant is to observe the effect 
of the different expansion history on dark matter structure formation [1]. In its standard version, 
quintessence is described by a minimally-coupled canonical field [2]. In this case scalar fluctuations 
propagate at the speed of light maintaining quintessence homogeneous even in the presence of dark 
matter clumps [3] . Quintessence can cluster only on scales larger than the horizon, where fluctuations 
have no time to propagate. However, observations on such large scales are strongly limited by cosmic 
variance and this effect is difficult to observe. 

A model of quintessence that can cluster on all observable scales has been recently proposed in 
[4, 5]. It is based on a single scalar degree of freedom with fluctuations characterized by a practically 
zero speed of sound. As explained in [4], there are several theoretical motivations to consider this 
case. In the limit of zero sound speed one recovers the Ghost Condensate theory [6], which is invariant 
under shift symmetry. Thus, there is no fine tuning in assuming that the speed of sound is very 
small: quintessence models with vanishing speed of sound should be thought of as deformations of 
this particular limit where shift symmetry is recovered [7, 8] . Moreover, using the tools developed in 
[7, 9], formulated in the context of an effective field theory, it has been shown that quintessence with 
an equation of state w < — 1 can be free from ghosts and gradient instabilities only if the speed of 
sound is very tiny, |c s | < 1CT 15 [4]. Stability can be guaranteed by the presence of higher derivative 
operators [7, 6], although their effect is absent on cosmologically relevant scales [4]. 

Apart from these theoretical considerations, a very important motivation to consider this model is 
that a series of galaxy and cosmic shear surveys are currently planned with the aim of understanding 
the nature of dark energy through its role in the structure formation. In this context the clustering 
scenario represents a phenomenologically interesting counterpart to the case of a smooth quintessence 
component. Indeed, quintessence with vanishing speed of sound actively participates to the formation 
of structures together with the dark matter and gives distinct modifications to the standard picture 
that can be strongly constrained by future data. 

In the past, several articles have investigated the observational consequences of a clustering 
quintessence in the linear regime, in particular, on the cosmic microwave background [10, 11, 12, 13, 14], 
galaxy redshift surveys [15], large neutral hydrogen surveys [16], the cross-correlation of the integrated 
Sachs- Wolfe effect in the cosmic microwave background with the large-scale structures [17, 18], or on 
weak lensing [19]. 

Theoretical investigations of the effect of dark energy on the nonlinear evolution of structures 
are particularly crucial. First of all, from linear theory alone it is difficult to distinguish the effects 
of quintessence on structure formation through its modification of the expansion history from those 
genuinely due to its perturbations. As we will see, the nonlinear evolution breaks this degeneracy. 
Furthermore, numerical simulations taking into account the gravitationally coupled evolution of dark 
matter particles and a clustering scalar field are still under construction and for the clustering scenario 
considered here they are totally missing. On the other hand, future redshift and weak lensing surveys 
will require very accurate predictions, both for the dark matter density and galaxy correlators, par- 
ticularly on nonlinear scales where the signal is larger. Finally, to conclude this series of motivations 
we remind that the study of nonlinearities is receiving a lot of attention in the context of primordial 
non-Gaussianities [20] . It is pertinent to ask whether a second clustering component could mimic the 
effect of primordial non-Gaussianities on the nonlinear evolution. 
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A first description of clustering quintessence in the nonlinear regime was given in [5] . There it was 
shown that in the limit of zero sound speed pressure gradients are negligible and, as long as the fluid 
approximation is valid, quintessence follows geodesies remaining comoving with the dark matter (see 
also [21] for a more recent model with identical phenomenology). In particular, reference [5] studied 
the effect of quintessence with vanishing sound speed on the structure formation in the nonlinear 
regime, in the context of the spherical collapse model (see [22] for a study of the spherical collapse 
when c 2 s of quintessence is small but finite). Due to the absence of pressure gradients, comoving regions 
behave as closed FRW universes and the spherical collapse can be solved exactly. The modifications 
to the critical threshold of collapse are small and the effects on the dark matter mass function are 
dominated by the modification on the linear dark matter growth function, which are also small. Today 
they are of the order of few per cent for realistic values of w. A larger effect occurs when one considers 
the total mass function, which includes the contribution of quintessence overdensities to the virialized 
halos. Indeed, quintessence contributes to the total halo mass by a fraction which increases at lower 
redshifts and is proportional to the ratio between quintessence and dark matter energy densities, 
i.e. ~ (1 + w) Q Q /Q m . 

In this paper we study the nonlinear regime of clustering quintessence using Eulerian Perturbation 
Theory (EPT). In particular, we extend the standard EPT approach for dark matter [23, 24, 25] to 
the presence of a second fluid, a clustering quintessence, comoving and coupled only gravitationally to 
dark matter. This is the first natural step to the study of nonlinear perturbations beyond the spherical 
approximation. In contrast to the spherical collapse model, this approach is perturbative and solutions 
can be found order by order. Notice that on small scales the EPT perturbative expansion for density 
correlators is not well defined, because it presents large cancellations between contributions of the 
same order. However, in the standard case it has been shown that classes of higher-order corrections 
can be resummed, leading to a well established perturbative scheme known, in its first formulation, 
as Renormalized Perturbation Theory (RPT) [26, 27, 28, 29, 30]. Complementary approaches can be 
found in [31, 32, 33, 34, 35]. 

In this work, we begin by considering the continuity equations for the dark matter and quintessence 
density contrasts and the Euler equation for their common velocity. Since gravitational observables 
are sensitive only to the sum of dark matter and quintessence fluctuations, we derive the continuity 
equation for the total density contrast. As both fluids are comoving, this equation and the Euler 
equation (together with the Poisson equation relating the total density to the gravitational potential) 
form a closed system in Fourier space, which can be solved perturbatively, as in the standard EPT 
approach. 

As in the standard case, the nonlinear couplings in the continuity and Euler equations are at most 
quadratic and the vertices are the same as those for a single dark matter fluid. The only difference 
is that the velocity diverg ence in the continuity equation is enhanced by the factor (1 + w^) £Yq / £l m . 
At linear order, this term is responsible for a rapid evolution of the growth rate at low redshifts, and 
changes the standard relation between the velocity divergence and the growth factor. Due to the 
absence of pressure gradients, the solutions for the linear growth and the linear growth rate of the 
total fluid can be written in integral form, as in the standard ACDM case. Using these solutions we 
are able to find simple fitting functions for these quantities, which generalize those currently employed 
in ACDM cosmologies [36, 37]. 

At higher order in the perturbative expansion clustering dark energy is responsible for an additional 
time-dependence of the kernels F n and G n defining the n-th order nonlinear corrections. The effect 
on Fi and G2 is of the order of the ratio between quintessence and total density perturbations, 
~ $pQ/{5p m + Spq), and gives distinctive signatures in the higher-order correlation functions such as 
the bispectrum. In particular, the reduced bispectrum, whose expression at leading order in EPT is 
independent of the linear power spectrum normalization, presents corrections only in the clustering 
case. Analogous corrections have been found in the halo mass function from the contribution of the 
quintessence mass to collapsed objects [5]. 

It is not the first time that EPT is generalized to the presence of several components. For instance, 
in [38] EPT has been applied to the problem of following the nonlinear evolution of baryon and cold 
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dark matter perturbations evolving from distinct initial conditions and in [39, 40, 41, 42, 43, 44] 
to the study of nonlinear perturbations in the presence of massive neutrinos. For modified gravity 
models it has been used in [45, 46, 47] to calculate the nonlinear power spectrum and, in particular, 
in [46, 47, ?] to compute the matter bispectrum. Higher-order observables, such as the normalized 
skewness S3 = (5 3 }/(S 2 ) 2 , have been also studied in [49, 50, 51] in the context of modified gravity 
models, where variations up to ~ 10% have been found. 

This paper is organized as follows. In section 2 we present the equations of motion describing 
the coupled evolution of matter and quintessence perturbations. In section 3 we solve the linearized 
equations for the density growth factor and the density growth rate and we discuss the solutions and 
fitting formulae. In section 4 we discuss the perturbative solutions in EPT. In particular, we derive the 
second-order solutions for the density and velocity fields, while in section 5 we derive the lowest-order 
observables: the density tree-level power spectrum and bispectrum. As a practical illustration of these 
results, in section 6 we compare the signal-to-noise expected for the effect of clustering and smooth 
quintessence with respect to the ACDM case in ideal measurements of the density large-scale power 
spectrum, bispectrum and reduced bispectrum in a box of 1 h~ 3 Gpc 3 at redshift z = 0.5. Finally, we 
present our conclusions in section 7. 

In addition, we present in appendix A.l a discussion on the analogy between the scalar field and the 
perfect fluid, with a derivation of the continuity, Euler and Poisson equations in the regime considered 
in this paper. In appendix A. 2 we derive evolution equations for the vertices in the spherical collapse 
approximation at all orders and in appendix A. 3 we discuss the redshift-space distortion effects in the 
clustering quintessence case. 



2 Equations of motion 

We consider a flat FRW background universe with metric ds 2 = a 2 (r)(— dr 2 + dx 2 ), where r is the 
conformal time, and a generic perfect fluid a with energy density p a , pressure p a , and peculiar velocity 
with respect to the Hubble flow v a . The continuity and Euler equations in an expanding background 
read 

d ^ + 3H( Pa + Pa ) + V- [(p a +p a )v a ] = 0, (1) 

<9t p a + p a V dr J 

where % = din a/dr is the conformal Hubble time and <3? is the gravitational potential satisfying the 
Poisson equation, 

V 2 $ = AirGa 2 ^ (5p a + 35p a ) . (3) 

a 

As explained in appendix A.l, these equations are valid only on scales much smaller than the Hubble 
radius. Furthermore, they assume non-relativistic fluid velocities, v <C c. For small density and 
pressure perturbations velocities remain small, independently of the speed of sound of the fluid. For 
the particular case of a fluid with zero speed of sound, such as dust or clustering quintessence, pressure 
gradients are suppressed and fluid velocities remain small even in the nonlinear regime. 

For small velocities one can neglect the time derivative of the pressure in front of the pressure 
gradient on the right hand side of eq. (2), v d T p <C Vp. In this regime the speed of sound is simply 
the ratio between pressure and energy density fluctuations, 



and eq. (2) reduces to 



cl,s = ^Pa/Spa , (4) 



^ + Hv a + (V a ■ V)v a = -^Zl^ _ W . (5) 
dT pa + Pa 
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For a fluid with vanishing speed of sound the first term on the right hand side of the Euler equa- 
tion vanishes and the fluid follows geodesies. This is the case for both dark matter and clustering 
quintessence. Thus, in their growing solution dark matter and quintessence are comoving [5] and we 
can take their velocities to be the same, i.e. v m = vq = v. 1 

Let us define the density contrast S a as 5 a = 5p a /p a , where p a is the background value of the 
energy density. In terms of this quantity, the continuity equation (1) for dark matter and quintessence 
becomes, respectively, 

^ + V-[(1 + <U^]=0, (6) 
-^--3wH8q + V-[(1 + w + 6q)v]=0, (7) 

where w = Pq/Pq is the equation of state of quintessence. In the limit cq jS = 0, the Euler equation 
for the common dark matter and quintessence velocity is the same, i.e. 

+Uv + (v ■ V)v = -V$ . (8) 

OT 

These are the equations describing dark matter and quintessence in the regime cq )S = 0. Note that 
the nonlinear couplings are only quadratic, as in the standard, pure dark matter case. 

In EPT it is useful to define the velocity divergence 9 = V-v. In fact, assuming v to be irrotational, 
which is a good approximation up to shell-crossing, the peculiar velocity v is completely described 
by its divergence 9 [25]. In Fourier space, the continuity equations for dark matter and quintessence 
become, respectively, 

95 mk 

—^T + % = -<x(Q1, 02) Qq^mm ' ( 9 ) 
- ZwU8 Q % + (1 + w)6 % = -a(q u q 2 ) 9^5 Q& , (10) 

where 

a(gi,<f 2 ) = H 2~ , (11) 

with qi = \qi\. Here and in the following, an integral f <fiq\(ftq25D{k — q\ —5*2) is implied over quadratic 
terms with repeated wavenumbers, such as on the right hand side of eqs. (9) and (10). There are few 
terms in eq. (10) that differ from eq. (9). The second term on the left hand side of eq. (10), absent 
when w = 0, comes from the fact that the energy of quintessence does not scale as the volume, while 
the factor 1 + w in front of the third term comes from the fact that the 3-momentum density of 
quintessence sourcing the variation of the density contrast is proportional to pq + Pq. The nonlinear 
couplings on the right hand side of eqs. (9) and (10) are both expressed in terms of a(q\, (7*2) and are 
the same for dark matter and quintessence. This is because pressure gradients are absent in both 
fluids. 

We can now use the Friedmann equation 3H 2 = 8nG(p m + Pq) and Sp m = 5pQ = to rewrite the 
Poisson equation (3) as 

v^ = ln 2 n m U n + 8 Q ^), (12) 



1 During matter dominance quintessence energy density is negligible and the dynamics is dominated by the 
gravitational potential wells of the dark matter. Assuming that quintessence velocity has not unreasonably 
"extreme" initial conditions, it will be rapidly driven to the dark matter velocity values by the force term in 
eq. (5), while velocity differences decay with the expansion. This situation is thus different from that described 
in [38] to study the gravitationally coupled evolution of baryons and dark matter, where the baryon fraction is 
always non-negligible with respect to the dark matter. 
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where Q a = fi Q (-r) = p a /{p m + Pq) for both matter and quintessence. Using this equation to replace 
the gravitational potential on the right hand side of eq. (8), we can rewrite the Euler equation in 
Fourier space as 



^ + h h + 1 n m n 2 (5 mS + 5 qS = - m&) % % 2 , (i3) 



«s) = ^-3-2 • ( 14 ) 



where 

(?i + g2) 2 gi • 92 
2 9i 2 9l 

To an observer measuring cosmological clustering through its gravitational effect, there is no 
distinction between 5p m and Spq: the observer will be only sensitive to the gravitational potential 
sourced by the total density perturbation bp = 5p m + SpQ, via the Poisson equation (3). We will thus 
concentrate on the description of the total density perturbation bp. A convenient way of describing 
the evolution of perturbations is to define the total density contrast 

6=^ = 6 m + 6 Q ^, (15) 

Pm ^'m 

which takes into account both dark matter and dark energy perturbations. With this definition the 
Poisson equation reads 

V 2 $ = 47rGa 2 p rn 5 . (16) 

An inadvertent observer could interpret S as due to dark matter perturbations only. However, as we 
will see below the time evolution of 5 when quintessence clusters is very different from the standard 
smooth case. Note that we could have defined 5 by dividing the total density perturbation Sp, not by 
the background matter energy density, but by the total background energy density p = p m + pQ. This 
would have had little consequence. Indeed, only the combination p m 5 is an observable via the Poisson 
equation so that the only effect, apart from a trivial rescaling in the evolution equations, would have 
been to change eq. (16) into V 2 <3? = AnGpS. Our definition has the advantage of recovering the usual 
expressions when quintessence perturbations are set to zero. 2,3 

By noticing that d r (£lQ/Q m ) = —3wHQ,Q/£l m , we can combine eqs. (9) and (10) into a single 
equation for S and write a closed set of equations involving just the total density contrast 8 and the 
velocity divergence 9, 

_1 + C9 % = -a(gi, £)%<% , (17) 
d9 3 



dr 



where we introduced 



C^^l + il + w)^ (19) 

As expected, the effect of quintessence perturbations on the growth of the total density contrast is 
proportional to (1 + w)Q,q/^L„^ so that it vanishes at early times, when VLq/Vtm — > or for w = — 1. 
Indeed, from eq. (13) the effect of 5q on the common velocity divergence is proportional to VlQ/Q, m . 
Furthermore, from eq. (10) the feedback of 6 on the density contrast is proportional to 1 + w. Note 
that the sign of the effect depends on the sign of 1 + w [11, 4]. As inside an overdensity 9 < 0, the 
total density contrast increases faster when 1 + w > 0, while increasing slower in the opposite case. 
For C = 1 we recover the standard case of canonical scalar field quintessence with c 2 s = 1, for which 
dark energy perturbations propagate as acoustic waves at the speed of light and quintessence remains 
smooth within a Hubble patch [1]. In this case there is no distinction between 5 and S m . 



2 While the matter density contrast has, by definition, the lower bound 5 m > — 1, for the total density contrast 
we have instead 5 > — [1 + Oq(t)/O to (t)]. 

3 To capture the change in the gravitational potential induced by dark energy perturbations, Ref. [52] defines 
a parameter Q such that the Poisson equation reads V 2( f> = 4-nGQp m 8 m . In our notation Q = 5/6 m . 
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Figure 1: The quantity C(z) as denned in eq. (19), as a function of redshift for several values of w. Unless 
otherwise specified, here and afterward we will denote w — —0.8, —0.9, —1.1 and —1.2 respectively by red 
dotted, red short-dashed, blue long-dashed and blue dotted-dashed lines. Continuos black lines will denote 
ACDM. The vertical line indicates the redshift of equality between matter and quintessence for a ACDM 
cosmology, i.e. z eq — 0.37 for the assumed fl m ,o — 0.279. 



Since the function C captures all the modifications to the equations of motion in the clustering 
case, it is useful to plot such quantity. In Fig. 1 we show C{z) for several values of w as a function 
of redshift. The vertical line indicates the redshift of equality between matter and the cosmological 
constant in a ACDM cosmology, that is z eq ~ 0.37 for the assumed f2 m> o = 0.279, where the index "0" 
denotes quantities evaluated today. 



3 Linear theory 

We will now study the linear solutions of eqs. (17) and (18) by neglecting the quadratic terms on the 
right hand sides of these equations. In Fourier space, we can write the linear solutions as 

Jf (r) s D(r)Sf , (20) 
" W/W fHf . (21) 



where D is the linear growth function. To derive the second equation we have used the linearization 
of eq. (17) and introduced the linear growth rate f as 

'=^- (22 » 

a In a 

a useful quantity to relate the velocity to the density perturbation in linear theory. Note that on the 
right hand side of eq. (21), unlike the usual case, the function C(z) is present at the denominator. 
Indeed, this comes from the factor in front of 6 in the continuity equation for the total density contrast, 
eq. (17). 

We will now study the evolution of D and / in the clustering quintessence case and we will 
derive fitting formulae for these two quantities, generalizing to the case of clustering quintessence 
the expressions employed in ACDM cosmology. For simplicity we will assume a constant w. See 
appendix A. 3 for an extension of the linear theory to redshift space. 
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3.1 Linear growth function 

By combining eq. (17) and the time derivative of eq. (18) at linear order one obtains the evolution 
equation of the linear growth function D. As a function of the scale factor a this reads 



d 2 D 



+ 



I, n , dlnC 

(i-3wn Q )- 



din a 2 

where the time derivative of C is given by 

dlnC 



din a 



dD 
din a 



- -n m CD = o , 



din a 



3w(C- 1) 
C 



(23) 



(24) 



For C = 1 we recover the evolution equation of the growth function in the case of a smooth quintessence, 
derived in [1]. 

Assuming that the scale factor today is do = 1, we can use the Friedmann equation to write the 
Hubble rate H = d In a/dt, where t is the cosmic time, as 



H(a) = H 



+ ^Q, a- 3(1+tu) 



1/2 



Then, deriving this equation with respect to the scale factor one obtains 

dlnH 3 
am a 2 

where Vt m = £l m $/{H 2 a?). Using the relation 

valid only for ACDM and for clustering quintessence, we can rewrite eq. (23) as 



d 2 D 
din a 2 



+ 



2 + 



dlnH dlnC 



din a din a 



dD dlnH n 

+ — D 



din a din a 



. 



(25) 



(26) 



(27) 



(28) 



We stress that this equation and the results that follow are valid only in the case of quintessence with 
zero speed of sound (or, in the limit where C = 1, for ACDM). Written in this form, it is easy to check 
that the evolution equation of D has two solutions. One is a decaying mode, D- oc H. The growing 
mode can be written in integral form as 



D + {a) 



H fi m ,o H(a) 



C{a) 



o [aH{d)f 



da , 



(29) 



where we have normalized D + in such a way that at early time, during matter domination, D + is 
equal to the scale factor a. This is the first of the main results of this work. Although derived for a 
constant w, one can check that eqs. (28) and (29) hold also when w depends on time. Alternatively, 
an explicit solution can be found in terms of hypergeometric functions by means of eq. (3.194) of [53]. 
This yields 



D + (a) 



y/l + x 



6w ' 



1 



-x 



where 



x(a) = 



6w ' 



5 — 6u> 

1 — ^m,0 



1 



6w 



,2 



5 

6w ' 



a 



-3w 



(30) 



(31) 



As already remarked, the solution (29) does not describe the case of a smooth dark energy com- 
ponent with w — 1. Indeed, such an integral solution relies on the absence of pressure gradients and 
on the fact that dark matter and quintessence move together along geodesies. In this case — and only 



8 



in this case — each comoving region evolves as an independent, unperturbed FRW universe. We can 
take advantage of this fact to re-derive eq. (29) in an alternative way, similarly to what is commonly 
done in the ACDM case [54] (see also [23, 55, 56]). Indeed, a spherical overdensity of dark matter and 
clustering quintessence of radius R can be described by the Friedmann equation for a closed universe 
[5], 

R\ 2 8irG , , K 



RJ = ^~^ m + ^)"^2' ( 32 ) 

where K is the curvature constant and the dot denotes the derivative with respect to the cosmic time. 
Let us describe deviations of the radius R from the scale factor of the background universe a by 

a = 1 - R/a . (33) 

Linearizing eq. (32) we obtain 

2Ha = -H 2 n m 5 + ^ ; (34) 

using the relation H = — (3/2)H 2 Q m C, which can be derived by taking the time derivative of eq. (25), 
we can rewrite this equation as 

(3dC) H = H5+ ~C . (35) 

_■ CI 

We can relate a to 5 by using the continuity equation inside the overdense region, i.e. [5] 

P + S^(P + Pq) = 0, (36) 

where p = p m + pq and the quintessence pressure Pq is unperturbed due to the absence of pressure 
gradients. Linearizing this equation one obtains 5 = 3 a C, which can be used to replace the parenthesis 
on the left hand side of eq. (35). With this replacement eq. (35) can be easily integrated to yield 

Ha) = \ H( a )K [^L-io, , (37) 

which presents the same time evolution as eq. (29). In this way we have consistently recovered the linear 
growing solution by linearizing the spherical collapse model, assuming that spherical overdensities 
behave as closed Friedmann universes. 

Now we turn our attention to the individual matter and quintessence perturbations. We introduce 
the linear growth functions of dark matter D m and quintessence Dq as 

These components do not evolve independently. However, once the solution for D is given, their 
evolution can be found by linearizing eqs. (6) and (7), to obtain 

dlna Cdlna' V ' 

dD Q l + w dD 

— 3wD Q = — — — — . (41 

a In a C a In a 

Notice that, from eqs. (40), (41) and (29), after few manipulations we can write D m {a) and Dq{cl) in 
integral form as a function of D(a) as 



D m {a) 

'o 



5 3D(a) 



2 2a 



£l m (a) da , (42) 
5 3D(a)~ 



2 2 



n Q {d)da. (43) 
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Figure 2: Upper left panel: ratio of the linear growth factor of the total perturbation, D + , in the clustering case 
to the smooth quintessence case. Upper right panel: the same as in the left panel for the matter growth factor, 
D m + . Lower left panel: ratio of the linear growth factor of the total perturbation D + in the clustering (thick 
line) and smooth (thin line) cases to the ACDM case, for fixed O m ,o today. Lower right panel: the same as in 
the left panel for D m +. 



In matter dominance, when D + = D m ^ + oc a, the growing mode for quintessence simply reads [4] 

D Qi+ = D ™,+ (matter dom.) , (44) 

so that 

D+ ( 1 + w fin \ . . . 

— = 1 + 7^- (matter dom. . (45 

a \ 1 — owilmj 

When quintessence dominates this is a poor approximation. Below, we will present a much better 
approximation, valid during both matter and dark energy domination. 

The effect of a clustering quintessence on the evolution of the total linear growing solution D + 
and of the linear dark matter growth function D m<+ is shown in Fig. 2, in the left and right panels, 
respectively. In the upper left panel we plot the ratio between D + in the clustering and smooth 
quintessence cases for several values of w. In the upper right panel we do the same for D m)+ . As 
expected, quintessence with zero sound speed enhances the clustering when w > — 1, while clustering is 
hindered for w < —1. As the effect is proportional to (l + w)0,Q/D, m , it increases in time and vanishes 
for w = —1. While the effect on the dark matter fluctuations is below the percent level for —1.2 < 
w < —0.8 (and for the assumed value £l m ,o = 0.279), the growth factor for the total perturbations 5 
is affected by modifications as large as 10% at redshift zero for models with |1 + w\ — 0.2. 
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In the lower left panel we plot the ratio of D + in the clustering (thick line) and smooth (thin line) 
cases to the ACDM case, for two values of w. In the lower right panel we do the same for D mt+ . 
In this case two opposite effects are into play, as one can clearly see from the bottom left panel of 
Fig. 2. On one hand, when w > — 1 quintessence dominates the Universe earlier than a cosmological 
constant, anticipating the accelerated expansion phase. This has the initial effect of suppressing the 
growth of fluctuations. However, since clustering quintessence contributes to the total perturbation 5, 
its effect on the linear growth function becomes important at low redshift, winning over the one of the 
accelerated expansion. As a consequence, the evolution of the gravitational potential is modified with 
respect to the smooth case and this distinctive signature can be used to constrain the dark energy 
parameter using the CMB [10, 11, 12, 13] or by cross-correlating the integrated Sachs- Wolfe effect with 
the large-scale structures [17, 18]. The effect of clustering quintessence on the matter linear growth 
is, on the other hand, much weaker, as shown in the bottom right panel. 

3.2 Linear growth rate 

From eq. (23) we can derive an evolution equation for the growth rate /, defined in eq. (22), given by 



din a 



1/ ~ x dlnC 

(l-3wn Q )-—- 
2 dma 



f = l^mC . (46) 



We remind the reader that the linear growth rate relates the linear velocity divergence 9 lm to the 
linear density contrast 5 lm as 

9 Vm = -H^ S Vm . (47) 

Thus, the significant quantity is given by the ratio f/C rather than / alone. Its explicit growing 
solution can be found using the definition (22) together with the integral solution for D + , eq. (29). 
This yields 



/+(«) 



Q m C , (48) 



where the ratio D + /a is given by eq. (30). The decaying solution is obtained from eq. (22) with 
L>_ oc H and from eq. (26), 

f_(a) = -hl m C. (49) 

In Fig. 3 the left panels show the effects of clustering quintessence on /+. In particular, in the 
upper left panel we plot the ratio of the growth rate f+ for clustering quintessence to the same quantity 
for a smooth quintessence component, for several values of w. In the lower left panel we plot the ratio 
of / + in the clustering and smooth cases to /+ in a ACDM cosmology. At high redshift the background 
evolution dominates the growth, as already discussed. At lower redshift /+ changes rapidly with time 
due to the clustering of quintessence. At z = the effect is quite large, of the order of ~ 20% for 
|1 + w\ = 0.1. 

In the right panels of the same figure, we show the effects of clustering quintessence on the ratio 
f+/C. As discussed in more details in appendix A. 3, it is this quantity, and not /+ alone, that describes 
the effect on redshift distortions at linear level. We show the ratio between f+/C for clustering 
quintessence and /+ for a smooth quintessence component (upper right panel) and for ACDM (lower 
right panel). As one can see, the growth of C in the ratio f+/C is very important, compensating — and 
marginally dominating — the growth of /+. Note that, from eq. (48), increasing the growth factor 
corresponds to decreasing the ratio f+/C. Indeed, quintessence with zero sound speed reduces f+/C 
for w > — 1, favoring the clustering, while the opposite happens when w < — 1. 

This qualitative behavior can be understood by introducing the growth index 7, so that /+ = QZi, 
where 7 is generally a time-dependent quantity. From eq. (46) and using 

= 3wn Q n m , (50) 

am a 
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Figure 3: Upper panels: ratio of the growth rate /+ (left) and of the quantity f+/C (right) for clustering 
quintessence to their counterparts in the smooth quintessence scenario (i.e., in both cases, /+) as a function 
of redshift for several value of w. Bottom panels: ratio of the / + (left) and of the quantity /+/C (right) for 
clustering (thick lines) and smooth quintessence (thin lines) to /+ in the ACDM case. 



the evolution equation for 7 in terms of the matter abundance Q m is given in our case by 

3wn Q inn rn ^- + \ + ni + 3w ( 7 - ^jn Q + 3w{c c ~ 1} - 3 -n^c = o . (51) 

Close to matter dominance, one can solve this equation as an expansion in Qq. In the smooth 
case, derived by solving eq. (51) with C = 1, one finds the growth index computed in [1], i.e. 

_ 3(1 -w) 3(l-w)(2-3w) , n , n 2s 
Tsmooth - -j^- + 2(5 _ _ 12u ,)"Q + ° W (52) 

~ 0.545 + 0.007 n Q + (1 + w) (0.025 + 0.005 fig) + . . . , 

where in the second line we have expanded as well with respect to 1 + uu. This can be compared to 
the growth index in the clustering case computed by solving eq. (51) with C given by eq. (19), i.e. 

6w 2 3w(72w 4 - 48w 3 - 3iw 2 + 56w - 25) 2 
7 ~ 5-6w (5 - 6u>) 2 (5 - 12w) Q+ ^ Q > (53) 

~ 0.545 + 0.007 Qq - (1 + w) (0.793 + 0.443 fig) + . . . . 
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Note that in the first line of eq. (52) we have corrected a typo in the denominator of eq. (Bll) of [1]. 
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Let us stress that 7 refers to the growth index of the total density contrast, which is different from 
the growth index of the density contrast of dark matter only computed, for instance, in the clustering 
case in [57] . 

The values of 7 sm ooth and 7 given in eqs. (52) and (53) are consistent with what shown in the upper 
left panel of Fig. 3, as one can check by approximating the behavior of the ratio between / + in the 
clustering and smooth cases with fi^ 7smooth , where 7 — 7 sm ooth = —(1 + io)(0.818 + 0.448 £Iq). Note 
also that 7— 7 smoo th strongly depends on Qq. Indeed, whereas in the smooth case the time dependence 
of the growth index is very weak, so that one can consistently approximate it as a constant also at 
low redshifts [58], in the clustering case this time dependence is more severe and such approximation 
is not viable. 

An alternative description of the growth function could be then given in terms of a reduced growth 
index 7 re( j, defined by f+/C = Qm ed , which presents a weaker dependence on cosmology and redshift. 
Indeed, one can check that 7 re d — 7smooth = (1 + to) (0.181 + 0.053 Oq), which agrees with what shown 
in the upper right panel of Fig. 3. An analogous behavior for the growth rate is discussed in [59] in 
the context of modified gravity models. 5 

Previous works (for instance [57, 14, 19]) have studied the clustering of quintessence concentrating 
on its effect on the matter perturbation S m , instead of the total perturbation 5. To connect to these 
works, we define the linear growth rate for matter perturbations as 



From eq. (40), it is related to / by 



_ d\nD m 

J™ = ~^P\ • ( 54 

a ma 



D m f m = D^. (55) 



3.3 Fitting functions 

As mentioned above, when quintessence dominates the Universe 7 becomes strongly time-dependent 
and the parameterization /+ = Qm, with 7 given in eq. (53), becomes a poor approximation (see, for 
instance, [60]). We have found the following approximation to be accurate in the range 0.1 < Q m < 1, 



U = c 



(56) 



As shown in the right panel of Fig. 4, for Q m > 0.2 the accuracy is better than 1% level for —1.15 < 
w < —0.85. For w = — 1 this formula reduces to the well-known fit given in [36] for ACDM. Following 
this reference, inverting eq. (48) and using in this equation the above approximation for /+, eq. (22), 
it is possible to find the fitting formula for the growth function 



a 2 



(57) 



As shown in the left panel of Fig. 4, for Q m > 0.2 and w > —1.5 the error is again less than 1%. When 
w = — 1 the commonly used fitting formula given in [37] for ACDM is recovered. 



5 Di Porto and Amendola study in [59] the growth of structures in a modified gravity model, parametrizing 
the growth rate /+ = ^7° A (1 + ?7da)- Their case is similar to ours. Indeed, in our case 7 rcc j and (1 + w)VlQ/£l m 
play the role of 7da and ?7da, respectively. Taking ??da to be constant, the authors of this reference constrain 
the values of 7da and tjua from a set of galaxy and Lyman-a observations, mostly at high redshift (z > 2). 
They find 7da = 0.6±g;| and ?7 DA = 0.0±g;| at l-a CL. In our model for |1 + w\ ~ 0.1, (1 + w)Vt Q /Vt m is 
very close to zero at z ~ 2, while it grows to ~ 0.2 at z = 0. Note, however, that due to the strong redshift 
dependence of (1 + w)f2g/f2 m it is difficult to extrapolate the validity of their results to our case. 
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4 Non-linear evolution 

The effects of a dominant quintessence component at low redshift are clearly not entirely captured by 
the linear growth factor. In particular, as quintessence clusters on all observable scales, we expect it 
to significantly affect the nonlinear evolution of structures. Thus, a proper description of the density 
correlators in the mildly nonlinear regime is required. Indeed, current and future redshift and weak 
lensing surveys will target this range of scales with great accuracy. 

As a first step in this direction, in this section we discuss the equations of motion (17) and (18) in 
the framework of EPT. In particular, we will derive their second-order solutions and use these to derive 
the tree-level expression for the bispectrum of the total density fluctuation 8, a valid approximation at 
large-scales. We leave for future work the study of nonlinear corrections to the density power spectrum 
and the bispectrum. 



4.1 Pert ur bat ive expansion 

In perturbation theory, it is useful to rewrite the fluid equations (17) and (18) using the growing 
function D + as time. In particular, we define r\ = \n.D + . Moreover, it is convenient to rescale 9 by 
defining 

e^-^e, (58) 

such that, according to eq. (47), at linear order & = 5, i.e. 

©f( V ) = ofirj) = D + (n)6f . (59) 

Note, again, that the definition of 0, eq. (58), contains also a factor C at the numerator, which is 
absent in the standard case. 

In terms of 0, the evolution equations (17) and (18) become 



dn k C 



j?2 V^fc °k> Q ^11^12 ' 

where we have used eq. (49) to replace — (3/2)il m C by /_. Written in this form, the continuity and 
Euler equations are the same as those derived for a smooth quintessence or ACDM cosmology [25] 



14 



except that the kernels a and (3 multiply the time-dependent function C~ l . Thus, the non-standard 
behavior of clustering quintessence is now encoded in the time dependence of the nonlinear couplings. 
The solutions to these equations can be written in perturbative form as 



oo 



§ = J2S^, Q = J2® {n) , (62) 
n=l n=l 

where 5^ = 5 lm and = lm . The n-th order solutions are proportional to n-th powers of the 
initial density 5 m and can be parameterized in terms of the kernels F n and G n and of powers of the 
linear growing mode solution D + as [25] 

4%) = F n (q lt q n ; V ) Z>» fa) 5%... 5% , (63) 
&f\ V ) = CM, . . . , q n ; r,) Dl(r,) 6% . . . 61 , (64) 

where F n and G n are homogeneous functions of degree zero of the momenta q\, . . . ,q n , with F\ = 
G± = 1. In these expressions a multiple integration over the momenta J d 3 qi . . . d?q n 5E>(k — YH=\ Qi) 
is implied on the right hand side. In the special case of matter dominance D + = a and F n and G n 
become time-independent and can be constructed from algebraic recursion relations [61]. In general, 
however, the kernels are time-dependent. 

As already pointed-out in section 2, since dark matter and quintessence are comoving, the system 
of equations involving 5 and (or 6) can be closed and we do not need to separately study the 
dark matter and quintessence to compute their evolution. However, it is interesting to compare the 
evolution of 5 with that of the dark matter and quintessence density contrasts, respectively S m and 
5q. Using the definition of 0, eq. (58), eqs. (9) and (10) can be rewritten as 

d5 m,k 1 ._ a(gi,g*) Q s _ (65) 

d5 Q,k 3u> 1 + ^ Q a(gi,g 2 ) Q x 

~W " U S Q* " ~C~ @ * = —c~ @ ^ Sq ^ ■ (66) 
One can then write the solutions to these equations perturbatively, similarly to what done above. 

4.2 Second-order solutions 

We now derive the second-order growing solutions in perturbation theory for the total density contrast 
5 and the velocity divergence 0. Thus, we replace on the right hand side of eqs. (60) and (61) the 



iluna t ?lin^ _ n f^Ain n „A aim/ 

over the momenta, 



linear growing solutions S^ n (rj) = D + {r])5"? and 0i n (r/) = D + (i])5^. This yields, after symmetrization 



d6 T (2) D\ ■ ■ 

- ©f - % - 4 2) ) = 4" 5 t . (68) 

where a s (qi,q 2 ) = [ct(qi,q 2 ) + (^(,02, 9i)]/2 is the symmetrized projection of a(qi,q 2 ). 

The solutions to these equations are usually parameterized in terms of F2 and G 2 , defined by 
eqs. (63) and (64) for n = 2 as 

4%) = F 2 (q u q 2 ; , (69) 

©f( V ) = G 2 (q 1 ,q 2 ;r ] )Dl( V )8 i i5l , (70) 
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so that the evolution equations (67) and (68) become 

01' 2 _„ „ 



dr] 
dG 2 
dr/ 



+ 2F 2 
+ G 2 - 



G 2 
/- 



C ' 
2 (G 2 - F 2 ) 



C 



One can verify that the solutions to these equations can be formally written as 



F 2 
G 2 



v 



v 



" r 3 «, + 2/3 D_(rj) ^_ v 2a s -2(3 
5 L>_(r/) 6 5 



e v-v 



3« s + 2/3 _ 3 f_(r,)f + (fj)D-(r,) ^ 2a s -2/3 
5 2f_(f,)f + ( V )D_(fi) e 5 



(71) 
(72) 

(73) 
(74) 



where we have used D + (rj) 



Note that these solutions take the same form as in the smooth case, 



except for the factor C(fj) at the denominator inside the time integrals, which in the smooth case is 
absent. 

Before discussing the clustering case, let us consider more carefully the smooth case. For C = 1 
the terms that do not contain the linear decaying solution D_ can be explicitly integrated to give a 
constant in time. Instead, the terms containing D- are generally time-dependent, so that each kernel — 
F 2 and G 2 — is characterized by a function expressing their time-dependence. For instance, one can 
use the angular average of F 2 and G 2} defining the functions v 2 and \i 2 respectively as v 2 = 2{F 2 ) and 
[i 2 = 2{G 2 ) [73]. These quantities also represent the second-order vertices of the nonlinear density 
contrast in the spherical collapse model (see appendix A. 2 for a generalization to higher order in 
perturbation theory). During matter dominance these functions become constant, i.e. v 2 = 34/21 
and \x 2 = 26/21. Indeed, in this case /+ = 1, /_ = —3/2 so that D-(rj) oc e~ 3r? / 2 , and the integrals 
can be solved to give constant F 2 and G 2 . As we will see, for a smooth quintessence or in a ACDM 
cosmology, /_ and /+ present a very weak dependence on time and the functions v 2 and ii 2 are very 
well approximated by their constant Einstein-de Sitter values. 

However, due to the presence of C(fj) in eqs. (69) and (70), in the case of clustering quintessence 
we need a third time-dependent function to fully characterize F 2 and G 2 . We choose to parameterize 
this additional time dependence with the function e, defined as 



e(r]) = / df/e^ 1 



c(v) 



(75) 



such that it does not involve the decaying mode D_ and vanishes when C = 1. With this definition 
F 2 and G 2 read 



F 2 (qi,q 2 ;r]) = - 
G 2 (qi,q 2 ;r)) = - 



1 " e(v) ~ 
1 " <rj) - 



3z^(77) 
2 

3^2(77) 



a s {qi,q 2 ) + 



1 - <v) - 

1 " e(v) ' 



2 

M2(?7) 



P(qi,&), (76) 

m,q2). (77) 



The evolution equations for v 2 and ii 2 can be found by averaging over the angles eqs. (71) and (72), 
using that (a s ) = 1 and (f3) = 1/3. They read 



dv 2 2 



d^ 2 
dr] 



+ M2 



J + 



v 2 ) 



2 

3C' 



(78) 
(79) 



with the initial conditions set by their solutions in matter dominance. An integral form of these 
functions can be given by taking the angular average of eqs. (73) and (74) and one can verify that 
they depend both on the linear growing and decaying solutions. 
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Figure 5: The function e defined in eq. (75) as a function of redshift for several values of w. The black continuous 
line corresponds to the zero value of e in ACDM and smooth quintessence cosmologies. 



It is possible to relate e to matter or quintessence fluctuations. Indeed, using D + = e 11 one can 
formally integrate eq. (40). This yields 



so that the second equality in eq. (75) can be rewritten as 

-D m ,+ _ Dq,+ 

e - 1_ ^7-^^7- (81) 

Thus, e is the ratio between the linear quintessence perturbations and the total perturbations, e = 
5pQ l /5p lm . In Fig. 5 we show the function e as a function of the redshift for several values of w. As 
quintessence becomes a non negligible component of the Universe, e becomes larger. Today it is of 
the order of ~ 0.05 for \w + 1| = 0.1. As we will see, this translates into similar corrections to the 
tree-level expression of the total bispectrum. 

In Fig. 6 we show the effect of quintessence on the vertices V2 (upper panels) and p,2 (lower panels) 
as a function of redshift and for different values of w. In particular, on the left panels we show the 
solutions in the smooth case, where we can notice the relatively small corrections induced by the 
different cosmological background to the Einstein-de-Sitter constant solutions. On the other hand, in 
the clustering case, shown in the right panels, the corrections are significantly larger, of the order of 
5% for |1 + w\ ~ 0.1 for both V2 and //2- These corrections can be estimated using eq. (82). 

As e explicitly characterizes the contribution of quintessence perturbations to the total one, it is 
useful to estimate the deviation of and (12 from their values in the smooth case in terms of this 
quantity. The exact explicit dependence cannot be written down in a simple form. However, we find 
that a very good approximation is given by 

V2 = ^2,smooth - (1.669 - 0.205 w)e , p 2 = /U2,smooth - (1-3 + 0.033 w) e , (82) 

where ^2,smooth and /^2,smooth are computed in the smooth case. As shown in Fig. 7, for Q m > 0.2, the 
error made using these approximations is less than about 0.1% for —1.2 < w < —0.8. 

Let us rewrite a s (qi, q 2 ) and /3(<fi, q 2 ) by making more explicit the dependence on the scalar product 
Qi ■ Q2, where qi = qi/qi. From their definitions, eqs. (11) and (14), one finds 

as{qUq2) = 1 + ~{q~ 2 + q - 1 ) ' (83) 
m,q 2 ) = ^(^ + q A + (q l -q 2 ?. (84) 

2 \q 2 qij 
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Figure 6: Upper panels: vertex for the pseudo-matter perturbations in the spherical collapse approximation 
in the smooth (left) and clustering (right) quintessence scenarios. Lower panels: same plots for the vertex 



Using these expressions we can reorganize the kernels F2 and G2 as multipolar expansions, in terms 
of their monopole, dipole and quadrupole contributions [62, 63] 

«(«,*) - 1 + a - ^ (§ + f) - \ (1 - « - 1) [1 - m ■ . (») 
- f + <> - ^ (! + !) - 1 - ' ~ f ) I 1 - 3 «> • «» 2 1 • ( 86 > 

As already mentioned, the first term corresponds to the second-order evolution in the spherical collapse 
dynamics. The dipole of the middle term is due to the nonlinear transformation from following 
mass elements in the Lagrangian description to describing the dynamics of the fluids in the Eulerian 
formalism. In the standard ACDM and smooth cases, this term is time-independent. In the clustering 
case quintessence contributes to the mass of a clustered object. As this mass is not conserved — the 
energy density does not scale as the volume — this nonlinear transformation receives a time-dependent 
correction proportional to e. The last term is due to the tidal gravitational field. Indeed, the dipole 
can be rewritten as 

[1 - 3(& • q 2 ) 2 ] %<% = -3 (q\q{ - \^)S^ 2 - ^)<% , (87) 
where (k^k^ — |5 ij )(5r is the Fourier representation of the tidal gravitational field (d l d^ — |5 ij V 2 )$. 
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Let us now consider the nonlinear correction to the density contrast of dark matter only, 5 m . We 
can plug the linear growing solutions (59) and (38) on the right hand side of eq. (65). This yields the 
second-order continuity equation for matter perturbations, 



85 



(2) 



'm,k 1 n (2) _ D+ Dm 



dr, C @ k = C as{quQ2) **** 

Making use of eq. (40) we can rewrite the partial derivative with respect to 77 = log D + as a derivative 
with respect to logZ) m ,+- Then, after multiplying this equation by D m /D, we can use eq. (81) to 
rewrite the continuity equation as 

1 - (1 - e) 4 2) = Dl . «,(&, q 2 ) SfS™ . (89) 



The velocity divergence on the left hand side is now weighted by a factor 1 — e = D m /D. Note also that 
the 1/C time dependence in the vertex present in the evolution of the total perturbation (eq. (67)) 
is now absent. As a consequence, we can describe the full time dependence of F m)2 in terms of its 
angular average, v m , 2 = 2(F m>2 )- 

The solution to this equation can be parameterized as 



ft^ = FmA^fr, r,)D 2 m>+ ( V )Sl5f 2 , (90) 



with 



If. 3 \ _ _ 3 /„ v m>2 



Fm,2(q\,q2) = ( 1 - 2^ m - 2 ) a ^^) + 2 [} 2~) PiQiifa) ■ ( 91 ) 



The evolution equations for z/ m ,2, obtained by taking the average of eq. (89), reads 

du m o 



d log D. n 



+ 2v m>2 - (1 - e)nz = 2 • (92) 



Note that if we rewrite F m2 in eq. (91) in terms of multipoles as in eq. (85), the e correction in the 
dipole term is absent. As already mentioned, this is a consequence of matter conservation and its 
relation to the transformation from Lagrangian to Eulerian space. More precisely, to the fact that the 
Jacobian of the coordinate transformation from Eulerian position x to Lagrangian position q is given 
by matter conservation, p m (l + 5 m )d 3 x = p m d 3 q, as J = (1 + <5. m (x)) _1 . 
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5 Power spectrum and bispectrum 

In this section we consider the leading-order contributions in perturbation theory to the power spec- 
trum and bispectrum of the total density contrast 5, which provide good approximations on large 
scales. The power spectrum P(k, n) is defined as 

{S^S^^Soik + k^Pik^) . (93) 

Assuming Gaussian initial conditions, it is possible to derive a perturbative expression for the nonlinear 
power spectrum in terms of the expansion in eq. (62). We have 

(^^) = (^ n 4 n ) + (4 2) 4' ) ) + ^ in 4' ) ) + --- ■ ( 94 ) 

where the first term on the right hand side defines the linear power spectrum, (<5| n 5 l lf) = 5o(k + 

k')P\m(k,ri), while the second and the third ones constitute one-loop corrections of the same order, 
i.e. fourth-order in the initial field 5 m . 

In terms of the power spectrum we have 

P(k, r,) = P lin (k, r,) + Pi_ioop(&, r)) + ■ ■ ■ , (95) 
where the linear contribution can be written as a function of the initial power spectrum, ((5 1 -? <5?) = 

k fc" 

5o(k + fc')Pin(fc), simply as 

Pi in (M) = D\(r,)P hl (k) , (96) 

and where Pi-i op(&> v) includes the two one-loop contributions mentioned above. Thus, at leading 
order in perturbation theory the effect of clustering quintessence is encoded in the linear growth 
function D + . On smaller scales one needs to consider nonlinear corrections. A consistent evaluation 
of the one-loop corrections for the power spectrum requires both the second- as well as the third-order 
solutions for the density field. We leave this for future work. 

The second-order solution 5^ 2 \ on the other hand, provides the tree-level expression for the density 
bispectrum, the leading contribution in EPT. The bispectrum is defined as 

(8^ ( V )S^ (77)^(7?)) = S D (h + k 2 + h)B(h, k 2 , h; n) . (97) 

Using Wick theorem, the tree-level solution for the bispectrum is given by 

B{k 1 ,k 2 ,k 3 ;r ] ) = 2F 2 (h, k 2 ; rj) Pi in (fci, v) Plmfa, v) + 2 cyclic . (98) 

The time evolution depends on the linear growth factor D + via the linear power spectrum (96) and 
on the functions e(rj) and v 2 (r\) in the kernel F 2 , eq. (76). For instance, for equilateral configurations, 
i.e. k± = k 2 = ks = k, 

B(k, k, k; rj) = 6F 2 (k, k; 77) Pi in (k, n) (99) 
= 3 {! Ufa) - \ [1 - £(17)]} D%{n) Pl{k) , (100) 

where in the second equality we have used eq. (96) to express the time dependence of P\m(k, i]). 

To highlight the shape-dependence of the bispectrum it is customary to introduce a reduced bis- 
pectrum Q(ki, k 2 , ks) defined as 

ntu u 1 ^- B(k!,k 2 ,k 3 ;ri) 

Q(h,k 2 ,k 3 ;ri) = — — — — — . 101 

P{ki,7]) P(k 2 ,7]) + 2 cyclic 

In the tree-level approximation, this is independent of the linear growth factor, as one can verify by 
using eqs. (95), (96) and (98) in (101), which yields 

no, u 1. \ 2 ^kfe ^ 2 ; rj) P in (A:i)Pin(fc 2 ) + 2 cyclic 

Q(ki,k 2 ,k 3 ;ri) = . 102 

Pin(^l)Pin(fc2) + 2 Cyclic 
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Figure 8: Effects of clustering quintessence on the tree- level bispectrum B {left panels) and reduced bispectrum 
Q (right panels) for the total density contrast 8, in the equilateral configuration, as a function of redshift. At 
tree-level Q is k- independent, while B depends on k, chosen here to be k — 0.1ft.Mpc -1 . In the upper right 
panel the constant value 4/7 for an Einstein-de Sitter cosmology is also shown as a dotted black line. The lower 
panels show the ratio with the corresponding B and Q for a ACDM cosmology. 



In an Einstein-de Sitter universe (matter dominance), as F2 is time-independent the reduced 
bispectrum is independent of redshift. Even in the ACDM and in the smooth quintessence scenarios, 
it shows a very mild time evolution. By contrast, as already explained, in the clustering case the time- 
evolution of the functions e and V2 is important and Q significantly departs from the standard value. 
Thus, the reduced bispectrum is sensitive to clustering quintessence. In particular, while the effect of 
clustering quintessence on the power spectrum and bispectrum is expected to be strongly degenerate 
with U m and 0%, this degeneracy disappears in the reduced bispectrum. Note that in redshift space 
the bispectrum depends on F2 but also on the kernel for the velocity divergence G2 which is affected 
through fj,2 by corrections of similar magnitude (see appendix A. 2). 

In the specific case of the equilateral configuration the reduced bispectrum becomes 

Q[k, k, k; r?) = 2 F 2 (k, k;r]) = l Mv) ~ \ [1 - efa)] , (103) 

which reduces to the constant value of 4/7 at early times, during matter domination. Note that the 
normalized skewness, S3 = (5 s ) /(5 2 ) 2 , is simply given by S3 = 3f2, so that the effect of clustering 
quintessence on the skewness is, using eq. (82), of the order ~ — 5e. 

In Fig. 8, in the left panels, we show the tree-level bispectrum in the equilateral configuration 
B(k,k,k; z) for k = 0.1/iMpc" 1 , as a function of redshift. In particular, the lower left panel shows 
the ratio between B(k, k, k; z) in the clustering and ACDM cases, for various values of w. Two effects 
are into play: the evolution of the linear growth function D + and that of the second-order kernel F2, 
through the functions V2 and e. The evolution is dominated by the fourth power of the linear growth 
function D + , see eq. (100). Indeed, D + is responsible for the turnaround at low redshift, as one can 
check in Fig. 2. We do not include a similar plot for the ratio between B(k, k, k; z) in the clustering 
and smooth cases, which is dominated by the linear evolution. The right panels of Fig. 8 shows instead 
the corresponding reduced bispectrum Q. The evolution is entirely due to the second-order kernel F2 
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Figure 9: Upper panels: the total reduced bispectrum Q(ki, k2, k^, z) for k\ = 0.05 h Mpc -1 and k2 — 2k\ as a 
function of the angle 9 between the two wavenumbers at z — (left) and z — 0.5 (right). Lower panels: ratio 
with respect to the ACDM value. 



and is more significant at low redshift, with corrections larger than 5% at z = for |1 + w\ = 0.1. 
Again, we do not include a similar plot for the ratio between Q(k, k, k; z) in the clustering and smooth 
cases, which shows a very similar behavior. 

Other configurations of particular interest are the collinear ones, where the triangle formed by the 
three momenta is flattened, such that the three sides satisfy the relation k% = k\ + k2- Thus, the scalar 
product between k\ and k2 is +1 while the other two are —1. In this case the dependence of F2 on V2 
drops off and the reduced bispectrum becomes 



Q(ki,k 2 ,h; V ) = [l-e(r?)j 



" (fci/fca + k 2 /h) P in (fci) P in (k 2 ) - 2 cyclic 
Pin(h) P m (k 2 ) + 2 cyclic 



(104) 



where the time-dependence is entirely factorized in the (1 — e) term. Thus, for collinear configurations 
the effects of clustering dark energy is described in terms of e(rf) alone, while no effect is expected in 
the smooth case. 

Let us study the shape dependence of the corrections due to clustering quintessence. In Fig. 9 we 
plot the reduced bispectrum Q(k\, k2, k%; z) with fixed k\ = 0.05 /iMpc -1 and k2 = 2ki, as a function 
of the angle 8 between k\ and /c2- At z = (left panels) the effect is of the order of 5% for |1 + w\ =0.1 
roughly for all the configurations considered and degrades to a few percent at z = 0.5 (right panels). 
One can notice that the effect of clustering quintessence is roughly independent of 0. In general, in 
the tree-level approximation the effect of clustering quintessence is similar on different scales, i.e. for 
different values of k\ and /c2- However, we expect higher-order corrections to become relevant on small 
scales and change this behavior. Notice that analogous corrections to the matter reduced bispectrum 
in the DGP model [64] are instead strongly dependent on the triangle shape, see for instance Fig. 10 
in [46] and Fig. 10 in [47] for a comparison with N-body simulations. 
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6 Signal-to-noise 



It is important to get an idea of the size of the corrections induced by clustering quintessence com- 
pared to cosmic variance. Since there is no way of distinguishing between matter and quintessence 
perturbations in any observation relying on gravitational effects, in the following we will only consider 
the power spectrum and bispectrum of the total perturbation 5. 

The effect of clustering quintessence on the power spectrum is encoded in the linear growth function 
D + , see eq. (96). As the power spectrum at fixed redshift is simply a function of one variable, i.e. the 
wavenumber k, the comparison of the corrections induced by quintessence with the expected statistical 
error is relatively straightforward. In the case of the bispectrum, however, a few plots of specific sets of 
triangular configurations can be misleading. In fact, the statistical significance of an observable, such 
as the galaxy bispectrum in redshift surveys, can only be assessed when all measurable configurations 
are accounted for. In this respect, it has been shown that the cumulative signal-to-noise expected for 
the galaxy bispectrum in surveys of the size of the SDSS main sample is comparable to the signal- 
to-noise for the power spectrum when all triangles down to mildly nonlinear scales are included in 
the analysis [65] . This holds true also when the full nonlinear covariance of the correlators and the 
survey geometry is considered. Moreover, a joint analysis of power spectrum and bispectrum can 
improve the constraints on cosmological parameters, particularly for those parameters responsible for 
the amplitude of perturbations and when the dark energy equation of state parameter w is allowed to 
vary [66]. 

Before presenting our conclusions, in this section we consider some simple estimates of the signal- 
to-noise ratio that we expect for the effect of clustering and smooth quintessence on the power spectrum 
and bispectrum in their leading-order, perturbative approximations. Such estimates do not relate to 
any specific large-scale structure observable. Thus, they should not be directly interpreted in terms 
of "measurable" departures from a ACDM cosmology. Rather, they provide a helpful comparison 
between the power spectrum and the bispectrum as large-scale structure probes of quintessence. More 
broadly, they aim to motivate further research more closely related to current and future observations 
in weak lensing or redshift surveys. 

For this purpose we consider a cosmological volume of 1 /i~ 3 Gpc 3 , significantly smaller than the 
volumes that will be probed by future surveys, such as Euclid, BOSS or LSST. We choose z = 0.5 as 
an intermediate redshift. However, since the time-evolution of the effect of clustering quintessence is 
very important, see Fig. 2, the redshift dependence could be used to increase our ability to detect this 
effect. Furthermore, we limit our calculation to scales k < 0.2/iMpc -1 , where we expect the linear 
and tree-level approximations to be sufficiently accurate for our purposes. Note that excluding small 
scales and neglecting higher-order corrections in P and B, which are relevant in the mildly nonlinear 
regime, is a conservative stance. Indeed, similarly to what we have discussed for F2 and G2, higher- 
order kernels are further affected by quintessence clustering. Thus, we expect higher-order corrections 
to P and B in the mildly nonlinear regime to be sensitive to clustering quintessence. 

The cumulative signal-to-noise for the power spectrum is simply defined here as 



The sum is intended over the measurable wavenumbers k from the fundamental frequency kf = 2n/L 
(defined by the assumed volume V = L 3 ) to the smallest scale k max in steps of kf. By Pcj(k) we 
denote the power spectrum for a (smooth or clustering) quintessence cosmology while Pa corresponds 
to the power spectrum in the ACDM case. The power spectrum variance AP 2 (k) is approximated by 
its leading, Gaussian contribution, 




(105) 



AP 2 



(k) 




P 2 (k). 



(106) 
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In the bispectrum case we define, mutatis mutandis, 




2 



[ff Q (fcl,fc 2 ,fc 3 ) ~ B A (h,k 2 ,k 3 )] 2 



E 



(107) 



fc, 



-max 



kf<ki, k2,k2 



where the sum runs over all triangular configurations defined by the wavenumbers k\ , k 2 and k 3 smaller 
than or equal to k mSLX . The bispectrum variance, again in the Gaussian approximation, is given by 



where we assume the reduced bispectrum variance to be dominated by the error on the bispectrum, 
so that AQ(k 1 ,k 2 ,k 3 ) ~ AB(k 1 ,k 2 ,k 3 )/[P(k 1 )P(k 2 ) + 2 cyclic]. 

Therefore, for these simple estimates we neglect the non-Gaussian contributions to the correlators 
variance, which we expect to be subdominant in the range of scales considered here, k < 0.2/iMpc -1 . 
Most importantly, we neglect the effect of covariance between different triangular configurations. As 
we have already stressed, these computations have purely illustrative purposes and do not assume any 
specific observable. Thus, we neglect any shot-noise contribution to the variances and any effect of a 
survey geometry, naturally expected, for example, in redshift surveys. 

In Fig. 10 we show the cumulative signal-to-noise for the effect of smooth (upper panels) and 
clustering (lower panels) quintessence with respect to a ACDM cosmology, on the total power spectrum 
(long-dashed lines), bispectrum (short-dashed lines) and reduced bispectrum (continuous lines). Note 
that w = —0.9 and w = —1.1 lead to very close bispectrum signals with the two corresponding curves 
almost coinciding. The reduced bispectrum is considered only in the clustering case since in the smooth 
case the corrections are negligible. In particular, on the left panels the signal-to-noise is shown as a 
function of the smallest scale considered, i.e. fc max for w = —0.9 (thick, red lines) and w = —1.1 (thin, 
blue lines). In the right panels, instead, it is shown as a function of w, assuming k mSLX = 0.2 h Mpc" 1 . 

The dependence on k max of the signal-to-noise for the power spectrum and the bispectrum is 
quite different. This has two reasons. On one hand the signal-to-noise for an individual triangular 
configuration of the bispectrum is smaller than that of the power spectrum for a given wavenumber 
of the same order. Thus, when only large-scale modes are considered the cumulative signal-to-noise 
in the bispectrum is smaller because only a few triangles can be measured on those scales. On the 
other hand, as we include smaller scales the number of available triangles grows more rapidly than the 
number of available wavenumbers. This is the case for the signal associated to the measurement of 
the correlators [65] as for the signal from the corrections induced by clustering dark energy, as studied 
here. One can check that the scale at which the two statistics display comparable signals depends on 
the volume, as expected. 

An interesting point, illustrated by the lines for the power spectrum and bispectrum on the right 
panels of Fig. 10, is that the dependence of the signal-to-noise on the equation of state w strongly 
depend on whether quintessence is clustering or not. For instance, in the particular example at 
hand, the ideal constraints on w for clustering quintessence can be a factor a few larger than in the 
smooth case. Indeed, the power spectrum and bispectrum are both mainly sensitive to the equation 
of state through the linear growth function D + , and Fig. 2 shows that, at z = 0.5, D + for clustering 
quintessence is very close to the one in ACDM, contrarily to D + in the smooth case. For the reduced 
bispectrum, which does not depend on D + , we find the opposite. Indeed, in the smooth case the 
dependence of Q on w is very mild while in the clustering case it becomes stronger, due to the 
modifications discussed in the previous section. 




(108) 



with sb = 6, 2 or 1 for equilateral, isosceles or scalene triangles respectively. 

Finally, we can define a cumulative signal-to-noise for the reduced bispectrum as 




(109) 



24 



■4* 

03 



10.0- 
5.0- 



Cumulative signal-to-noise 
smooth Q vs. ACDM 
V= 1 /r'Gpc 3 

z = 0.5 ^ 



1.0 



0.5 U 



/, 

I / 

/ // 
»/ 
a 
a 



<', 



Power spectrum, P: 

— — w = - 0.9 

— — w = - 1.1 
Bispectrum, B: 

W = - 0.9 

w = - 1.1 



0.05 



0.10 

ax [h Mpc" 1 ] 



0.15 



0.20 



■4* 

V 



10.0- 
5.0- 



Cumulative signal-to-noise 
clustering Q vs. ACDM 
V = 1 ft-'Gpc 3 
z = 0.5 



1.0- 
0.5- 




Reduced bisp., Q: 

w = - 0.9 

w = - 1.1 



0.05 



0.10 

ax [h Mpc" 1 ] 



0.15 



0.20 



5.0 



03 



0.5 



^ i I I 1 1 1 1 I 1 1 I 1 

'Cumulative signal-to-noise: 

smooth Q vs ACDM 
■s \ 

\ 

\\ 


✓ 1 1 

// 
>'/ 


\ 


\ 1 






if 

'l 






z = 0.5 


— — p 




V = 1 ft, Lrpc 


B 


\ 
L 


^max= 0.2/lMpC" 1 


-1.2 -1.1 


-1.0 
w 


-0.9 -0.8 



5.0 



03 



0.5 



Cumulative signal-to- 


noise: 




clustering Q vs ACDM 






s 

vs. * 
\ \ 


✓ 

✓ 

/ 

/ > 
/ X 
/ / 
/ / 
/ / 

'/ ' 




V 


if/ 




- __, X 


1 






if v 

i !l 


z = 0.5 


B \ 

Q 


= 1 /r'Gpc 3 " 


r I £ - 

1 / nnax — 


0.2 hMpc- 1 



-1.2 



-1.1 



-1.0 



-0.9 



-0.8 



Figure 10: Cumulative signal-to-noise for the effect of smooth (upper panels) and clustering (lower panels) 
quintessence on the total power spectrum (long-dashed lines), bispectrum (short-dashed lines) and reduced 
bispectrum (continuous lines) for a volume of 1 ft, -3 Gpc 3 at z = 0.5. On the left panel the signal-to-noise is 
shown as a function of the smallest scale considered, i.e. /c max for w = —0.8 (thick, red lines) and w = —1.1 
(thin, blue lines). In the right panels the signal-to-noise for power spectrum, bispectrum and reduced bispectrum 
is shown as a function of w, assuming fc max = 0.2 h Mpc -1 . 



In Fig. 11 we show the signal to noise for the corrections induced by clustering dark energy with 
respect to the smooth case. For the power spectrum the expression is given by 

S_\ 2 k f? [P Qc (k) - P Qs (k)f mm 

where Pq c and Pq s represent the power spectrum for clustering and smooth quintessence, respectively. 
Analogous expressions can be written for the bispectrum and the reduced bispectrum. Note that the 
fact that the signal-to- noise for the power spectrum and the bispectrum in the right panel of Fig. 11 
are so close is a coincidence due to the choice of fc maX5 as can be seen from the left panel. 

The results of Figs. 10 and 11, particularly the relative size of the signal-to- noise of the power 
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Figure 11: Cumulative signal-to- noise for the effect of clustering quintessence on the total power spectrum (long- 
dashed lines), bispectrum (short-dashed lines) and reduced bispectrum (continuous lines) with respect to the 
smooth quintessence scenario assuming a volume of 1 h~ 3 Gpc 3 at z = 0.5. On the left panel the signal-to-noise 
is shown as a function of the smallest scale considered, i.e. k max for w = —0.9 (thick, red lines) and w — —1.1 
(thin, blue lines). In the right panel the signal-to-noise for power spectrum, bispectrum and reduced bispectrum 
is shown as a function of w, assuming fc max = 0.2 h Mpc -1 . 



spectrum versus that of the bispectrum, are strongly dependent on the volume. For instance, larger 
volumes favor bispectrum measurements as the number of configuration increases significantly. Fur- 
thermore, these results strongly depend also on the redshift through D + and e (see Fig. 2 and 5). 
Thus, one should be particularly careful in extrapolating them to other redshifts. For instance, as can 
be seen in Fig. 8, the choice of z = 0.5 corresponds to the specific case where the effect of clustering 
quintessence is maximum on the bispectrum, while it gets significantly larger at smaller z for the 
reduced bispectrum. 

For these reasons, what shown here does not provide an exhaustive comparison between different 
statistics or a complete picture of the signal expected. However, this analysis indicates that both 
correlators should in principle be equally sensitive to corrections induced by quintessence, both in the 
clustering and in the smooth case. 

The relatively high signal-to-noise expected in the ideal set-up assumed here motivates further 
studies. Indeed, an actual detection of such features would require the ability of distinguishing them 
from other nonlinear effects or, in the case of redshift surveys, from galaxy bias and redshift distortions. 
The detailed study of the detectability of possible departures from a ACDM cosmology in actual 
observation is thus left to future work. 

7 Conclusions 

In this work we considered the case of a quintessence characterized by a vanishing speed of sound. In 
this case quintessence perturbations grow on all observables scales, inducing relevant effects on the 
evolution of structures when the dark energy component comes to dominate the energy density of the 
Universe. 

At late time, both dark matter and quintessence perturbations act as a source for the gravitational 
potential and they are practically indistinguishable by gravitational observations. Based on this, we 
introduce a total density perturbation as a weighted sum of matter and quintessence perturbations. 
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Since quintessence is comoving with dark matter, the evolution of such a quantity is determined by 
a closed set of equations: the continuity equation for the total density field, the Euler equation for 
the common velocity field of the two components, and the Poisson equation relating the gravitational 
potential to the total density perturbation. This allows us to study the evolution of the total density 
fluctuation in Eulerian Perturbation Theory, in complete analogy with the usual treatment of matter 
fluctuations in a ACDM or smooth quintessence cosmology. 

The equations of motion for the total perturbations are equivalent to those for the matter per- 
turbation alone, with a simple correction: the linear term in the velocity divergence of the continuity 
equation is proportional to the function C = 1 + (1 + w)0,Q/U m . At early times, when quintessence 
is negligible, C = 1 and we recover the standard evolution. At late times clustering quintessence 
increases the growth rate of fluctuations. 

At linear order it is possible to obtain an exact integral expression for the growth of the total 
perturbation, eq. (29). This solution relies on the fact that for a quintessence with vanishing speed 
of sound comoving regions behave as independent FRW universes. This integral expression allowed 
us to derive simple fitting functions for both the growth function and the growth rate of the total 
fluctuations. 

Beyond linear order, our set-up allows to straightforwardly apply to the clustering quintessence 
scenario standard EPT, but also more efficient resummation techniques such as Renormalized Pertur- 
bation Theory [26, 27, 28, 29, 30] and the Renormalization Group approach [31, 32, 33, 34, 35]. 

We showed that linear theory does not fully describe the rich phenomenology of a quintessence 
with zero speed of sound. Indeed, we found significant effects on the late-time evolution of higher- 
order perturbations. These can affect the total power spectrum over a wide range of observable scales, 
where the evolution of perturbations becomes nonlinear. 

Since they directly depend on nonlinear corrections over the linear density field, also higher-order 
correlation functions, such as the bispectrum, are affected by the clustering of quintessence at low 
redshift. In particular, we studied second-order solutions in EPT for the density contrast and velocity 
fields. In terms of these solutions we derived the leading-order (or tree-level) contribution in EPT 
to the total bispectrum. On large scales, this is expected to be a good approximation to the fully 
nonlinear bispectrum. 

In particular, we showed that the reduced bispectrum, which is normalized in such a way as to be 
independent of the linear evolution, receives significant corrections only in the clustering case. These 
corrections are of the order of 5pQ/(5p m + ^Pq), i-e- the ratio between the quintessence and total 
density perturbations, which at z = amounts to 5% for |1 + w\ = 0.1. These signatures offer a 
practical way of distinguishing the clustering scenario from the smooth one with the next generation 
of redshift and weak lensing surveys such as BOSS or Euclid. Corrections of the same magnitude to 
the reduced matter bispectrum are expected as well from non-Gaussian initial conditions. However, 
notice that, at least on large scales, such corrections present a different (in fact, opposite) redshift 
evolution as well as different dependences on scales and shapes (see, for instance, [78]). 

In Section 6 we provided a simple estimate of the signal-to-noise ratio expected for the effect of 
quintessence on the power spectrum and bispectrum of the total density field. We limited our analysis 
to an ideal box of 1 /i~ 3 Gpc 3 at fixed redshift z = 0.5. In particular, for the linear power spectrum 
we considered the signal expected for the difference between the predictions of smooth quintessence 
and ACDM, Pq s — P\, of clustering quintessence and ACDM, Pq c — P\, and of clustering and smooth 
quintessence, Pq c — Pq 3 . We performed the same analysis for the tree- level bispectrum and the 
reduced bispectrum. In order to provide a fair comparison between power spectrum and bispectrum, 
we included all measurable triangular configurations down to a given k max . Below a given scale, the 
signal-to-noise ratio of the bispectrum becomes more important than the one of the power spectrum. 
Interestingly, this takes place at a larger scale for the clustering case. 

This preliminary analysis is clearly very limited. We discussed simply the density perturbations 
without considering a particular observable. Moreover, we did not include higher-order nonlinear cor- 
rections. These are expected to be relevant for wavenumbers close to the maximum value considered 
here, i.e. /c max = 0.2/iMpc -1 . Nevertheless, our results emphasize the importance of a joint analysis 
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of the power spectrum and bispectrum in future redshift and weak lensing surveys, possibly extend- 
ing over a large redshift range to take full advantage of the rich time-dependence in the clustering 
quintessence scenario. Furthermore, they suggest that smaller scales, where the evolution of perturba- 
tions is nonlinear, are likely to be affected significantly by quintessence clustering. Thus, extending this 
analysis to the nonlinear regime — somehow a necessary task — can significantly improve the constraints 
on dark energy with respect to those forecasted assuming only linear theory predictions. 

Given the significantly different phenomenology between the smooth and clustering cases, it is 
crucial to develop accurate theoretical predictions to compare with cosmological observations. As 
shown, the inclusion of higher-order perturbations, either as corrections to the power spectrum or in 
higher-order statistics, can be very important in this process. In this work we presented a first step 
in this direction using perturbative techniques based on Eulerian Perturbation Theory. We leave to 
future study the extension to more efficient resummation schemes. 
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A Appendices 

A.l Scalar field and fluid equations with pressure 

In this section we show that a scalar field with zero speed of sound of fluctuations satisfies eqs. (1), 
(2) and (3). 

Let us consider the action for fc-essence, i.e. [69, 70] 

S = J d A x^~g P(<f>, X), X = -sTd^A ■ (HI) 

The evolution equation of <j> derived from this action is 

V„( 2P X &*4>) + = , (112) 

where Pj = dP/df and is the covariant derivative with respect to the metric g^ v . The energy- 
momentum tensor of this field can be derived using 

T -_J^ll_ (113 ) 

and can be written in the perfect fluid form as [71] 

= (p + p)u^u v + pg^ v , (114) 

once we identify the rest-frame energy density and pressure and the 4-velocity of the fluid, respectively 
as 

p = 2P x X-P, p = P, „ M = -M. (115) 



Now we will show that the scalar field satisfies relativistic fluid equations. The relativistic conti- 
nuity equation follows from the evolution equation (112). Using the third equality in (115) to replace 
<9 M by the 4-velocity and multiplying it by y/X this equation reads 

-u^( 2P tX yfX)yfX + P^yfX - 2P X XVX = ■ ( 116 ) 
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By employing that y X = u^d^cf) this equation can be rewritten as 

[d^ 2P X X) - P x d^X - P+d^] + 2P X XVX = , (117) 

where, using the first equality in eq. (115), we recognize in the quantity in brackets the gradient of p, 
i.e. d^p. Thus, eq. (117) can be rewritten as a continuity equation of the energy, 

u^p+(p + p)V^ = 0, (118) 

which, alternatively, can be derived from projecting the conservation equation of the energy- momentum 
tensor along the fluid 4-velocity, u v XI ^ v = [72]. 

As shown in [5], the relativistic Euler equation is an identity for the scalar field (p. Taking the 
derivative of the definition of X in eq. (Ill), 

d v {d^d^) = -d v x , (ii9) 

and rewriting it in terms of the 4-velocity = —d^(f>/Vx we have 

2u< i VxV li {VXu v ) = -d u X , (120) 

and thus 

2Xu^V^u u = -(g VfM + //J^I . (121) 

In the parenthesis we recognize the projector on hypersurfaces orthogonal to u^. Then we can use 
d^P = P^d^cf) + P : xdfj,X to replace d^X in this equation, noting that the term proportional to d^cf) 
vanishes when multiplied by the projector orthogonal to u^. Thus, we obtain 

(p + pKVX = + u v u^p , (122) 

which, alternatively, can be derived from the conservation equation (g vp + u u u p )X7 ^T^ p = [72]. In 
conclusion, as expected the scalar field satisfies the relativistic continuity and Euler equations, and 
thus it is dynamically equivalent to a perfect fluid. As the speed of sound of this fluid is defined as 
c s = P,x/p,x [71], the right hand side of the Euler equation (122) can be rewritten as 

(p + p)u»V,u» = -c 2 s (g^ + u v vP)d» P , (123) 

which shows that in the limit c s = the fluid satisfies geodesic motion u p V IM u u = [5]. 

Let us neglect for a moment metric perturbations, which will be reintroduced later on. Consider 
a coordinate system (t, x l ) where the 4-velocity can be written as 

«" = 7 (1,« < ), (124) 

with v i = dx i jdt being the 3- velocity of the fluid and 7 = 1/ y/l — v 2 the relativistic factor. Note 
that the 3-velocity is related to the field by v = V<p/(j). Multiplying eq. (118) by 7 and using in this 
equation the definition (124), summing it to the v = component of eq. (122), one finds 

dt [l 2 {p + v 2 p)] + V • [ 7 2 (p + p)v] = , (125) 

where we recognize in r y 2 (p + v 2 p) the energy density in the rest frame defined by (t, x % ). Alternatively, 
this equation can be derived more directly from the v = component of the conservation equation of 
the energy-momentum, V^T^q = 0. 

Multiplying the v = component of eq. (122) by v and subtracting it from the spatial component 
of the same equation one finds 

l 2 {p + P) dtv + (v-V)v =—Vp — vdtp, (126) 

which can be also simply derived from V^T^* = 0, after using eq. (125). 
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Let us consider fluctuations around an equilibrium state of the fluid, characterized by time depen- 
dent energy density and pressure p(t) and pit) = wp(t) and velocity v = 0. Equations (125) and (126) 
simplify in the limit of small velocity, i.e. for v <C c = 1. In this case they can be written as 

d tP + V-[(p + p)v\=0, (127) 
(p + p) d t v + (v-V)v = -Vp-vdtp . (128) 

We can now reintroduce gravity in these equations by assuming a perturbed flat Friedmann metric, 

ds 2 = -dt 2 {\ + 2$) + a 2 (t)(l - 2$>)dx 2 , (129) 

with <E>, ^ <C 1. Gravity enters through the covariant derivatives of the 4- velocity on the left hand sides 
of eqs. (118) and (122). In the continuity equation it yields the term {p + p)u^V ^ In \f—g — (p+p)3H 
and in the Euler equation it introduces the term (p + p^^Y 1 u p ~ (p + p)(2Hu l + 9*$), where we 
have neglected subleading terms in v/c and Thus, with these new terms and using comoving 

coordinates with v l ~ au l , eqs. (127) and (128) become 

d t p + 3H{p + p) + ^V ■ [{p + p)v} =0, (130) 

d t v + Hv+ l (v- V)v + -V$ = — f -Vp + tf d t p | , (131) 

a a p + p \a J 

which are equivalent to eqs. (1) and (2) after introducing the conformal time r related to t by dt = adr. 

The gravitational potential $ can be related to the energy density and pressure perturbations 
by the Einstein equation. Neglecting time-variations of <I> and ^ of order H, which are small in the 
sub-Hubble scale dynamics [5], the 00 component of the Einstein equation yields 

V 2 ^ = AitGa 2 5p , (132) 

while the traceless part of the ij component yields 

V 2 ($ - = YlnGa 2 bp . (133) 
By combining these equations one obtains the Poisson equation (3), 

V 2 $ = 4irGa 2 (5p + 35p) . (134) 

Note that in deriving eqs. (130) and (131) we did not assume that energy and pressure perturba- 
tions are small with respect to their background value. In this sense, these equations are nonlinear 
in the energy density and pressure. However, for large pressure gradients one expects that the veloc- 
ity, sourced by the right hand side of eq. (131), becomes close to relativistic values invalidating the 
assumption v <C c. This does not happen when pressure gradients are suppressed by the smallness of 
the speed of sound, as for dust or clustering quintessence. Indeed, in this case the right hand side of 
eq. (131) vanishes, and these equations consistently describe the nonlinear regime. 



A. 2 Vertices in the spherical collapse approximation 

In the spherical collapse model one assumes spherical symmetry around x = 0. As a consequence, the 
linear density field 5\\ a {k) depends only on the norm of k, k = \k\ and the nonlinear solutions for the 
density contrast 5 and velocity divergence © can be written as [73], 

MT) = £^fWK\ (135) 

n=l 

oo , , 

e sc(r) = ^MI2 J D«( T)e «, (136) 
n. 

71=1 



30 



where v n and [i n are the angular averages of the kernels, 

v n = n!(F n (gi, ...,q n ))=n\ J ^ . . . ^F n (<fi, . . . ,q n ) , (137) 

Hn = n\{G n {qi, ...,q n )) = n\ J ^ . . . ^G„(<fi, . . . ,q n ) , (138) 

e = J d 3 k5i n (k) and i^i = /ii = 1. 

From the equations of motion (60) and (61) we find the recursive equations for the evolution of v n 
and n n , 

n-1 



dUr. 



du n 1 ^-^ I ' n\ 

— + n V n - jl n = — I I VmVn-m , (139) 
m=l ^ ' 

/ \ 3 f^^C . . 1 ^ — > / 71 \ . . 

+ (n- \)n n + ^—p2-{^n ~ Vn) = ^ 2^ I m J^mMn-m • (140) 

/+ m=1 V / 



3?? 

In particular, for n = 2 we have 

w +^-» = c- < 141 » 

3 m C 2 
■^-+^+2^-^-^ = 3C- (142) 

A. 3 Redshift distorsions 

The derivation of the Kaiser formula [74] can be straightforwardly extended to the case of clustering 
quintessence. We follow the presentation of [25] but we focus our considerations to the fluctuations 
in the number density of galaxies 5 g , rather than the matter overdensity 5 m . As both these densities 
are conserved, there is in fact no essential difference. In this way we directly refer to a generic biased 
population like the galaxy distribution. For simplicity, we work in the parallel-plane approximation. 
The mapping between redshift and position space, in comoving coordinates, is given by 

s = x + ^z, (143) 

where z is the direction along the line-of-sight. From the conservation of the number of galaxies we 
can find the relation 

(1 + 5 g , s ) d 3 s = (1 + 5 g ) d 3 x (144) 

between the galaxy overdensity in redshift space 5 9:S and the same quantity in position space S g . 
As the volume elements are related by d 3 s = J(x)d 3 x, where J is the Jacobian of the coordinate 
transformation from s to x, one finds 

m- ^y . ("si 

The Jacobian is explicitly given by J = |1 + V Z v z /T-L\. In Fourier space, under the assumption of 
V ' Z v z /T-L <C 1, this yields 



= / ^e- A ^ k ^l n [5 g (x)-V z v z (x)/H] . 



(146) 



This equation can be written in terms of © = —Cj (%/) V • v defined in eq. (58) (here for simplicity 
we use the convention f = f + ) as 

«*) = E u— tv r£- / d3qi • • • d3qMk ~ £ *) 

n=l l " 1J ' VO/ 3 i=l (147) 

— ©(<?2 (?n) , 

<?2 Qn 



C 
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where fi = k z /k and fa = q z ,i/qi- For C = 1 we recover the standard expression of [75]. At linear 
order this expression reduces to 



(148) 



Assuming a linear bias relation between the galaxy density and the total linear perturbations, 
i.e. Sg 111 = b\ <5 lm , and replacing lin using eq. (59), we can write, in the linear approximation 



din/ 



/ 



.2 din 



c 



f 



rlin / 



(149) 



This is the analog of the Kaiser formula, except that / is here replace by f/C. Using this expression, 
the linear galaxy power spectrum in redshift space is thus given by 

P g , s (k) = P g (k) {1 + (3 fi 2 ) 2 , (150) 

where we have introduced the parameter which in the case of clustering quintessence is given by 

/ 



(151) 



Thus, at linear order the corrections to the redshift distortions formula due to clustering quintessence 
enter only through the function C in the parameter (3. However, we can also express f3 in terms of the 
linear growth rate f m defined in eq. (54), and the linear bias parameter b m ,i, assuming a linear bias 
relation between the galaxy density and the matter linear perturbations, i.e. 5 g m = 6 m ,i ^m 1 - Using 
eq. (55), we can rewrite (3 in eq. (151) as 



fr, 



fm Dm 

D hi b m A ' 



(152) 



recovering the relation employed in [19]. 

From eq. (147) one can also derive an expression for the galaxy bispectrum in redshift space. 
Assuming a nonlinear local relation between the galaxy overdensity 5 g and the total density contrast 
5 and Taylor expanding it in terms of powers of S one obtains [76] 



(153) 



n=l 



Furthermore, the left hand side of eq. (147) can be expanded in terms of redshift-space kernels Z n as 
[77, 75] 

oo „ n 

5 g , s (k, V ) = J2 DT l(v) J d 3 q 1 ---d\ n 5 D (k-Y,^)Zn(qi,...,qn;v)Sf 1 ---Sl. (154) 



n=l 



Fourier transforming eq. (153) and using it to replace 6 g (k) on the right hand side of eq. (147) one 
obtains, up to second order, 

Z 1 (fc,7?) = 6 1 (l + /?/i 2 ), (155) 



and 



Z2(ki,k 2 ;rj) = b\ F 2 (ki, k 2 ; rf) + j3fi 2 G 2 (ki, k 2 ; v) 



Ml 



(i + ^) + £ (l + ptf) 



k 2 



(156) 



+ 



2 ' 



with F 2 and G 2 defined in eqs. (69) and (70). Thus, as before these expressions for Z\ and Z 2 can be 
quickly recovered from the standard ones after replacing / by f/C. 
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From eq. (154), the redshift-space, tree- level, galaxy bispectrum is given by 

B g>s (ki,k2,k3;ri) = 2Z 2 (ki,k2;rf)Zi(ki,r])Zi(k2,r])P^ n (ki,r))P^ n (k2,r]) + 2 cyclic. (157) 

The effects of clustering dark energy enters in the ratio f/C but also in the corrections to the second- 
order kernels F 2 and G 2 , which are both of the same order. In addition, one should carefully consider 
the linear and quadratic bias parameters b\ and b 2 defined in terms of the total density. One can 
expect for instance that for an ideal population of galaxies characterized by a conserved comoving 
number density, the corresponding linear bias would present a significant dependence on redshift at 
late times. The implications of the clustering quintessence scenario for galaxy bias are beyond the 
scope of this work and will be considered elsewhere. 
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